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Abstract. We study dynamical consequences of the kink instability of a twisted coronal flux rope, using the force-free coronal 
loop model by Titov & Demoulin ( 1999) as the initial condition in ideal-MHD simulations. When a critical value of the twist is 
exceeded, the long- wavelength (m =1) kink mode develops. Analogous to the well-known cylindrical approximation, a helical 
current sheet is then formed at the interface with the surrounding medium. In contrast to the cylindrical case, upward-kinking 
loops form a second, vertical current sheet below the loop apex at the position of the hyperbolic flux tube (generalized X line) 
in the model. The current density is steepened in both sheets and eventually exceeds the current density in the loop (although 
the kink perturbation starts to saturate in our simulations without leading to a global eruption). The projection of the field lines 
that pass through the vertical current sheet shows an S shape whose sense agrees with the typical sense of transient sigmoidal 
(forward or reverse S-shaped) structures that brighten in soft X rays prior to coronal eruptions. The upward-kinked loop has the 
opposite S shape, leading to the conclusion that such sigmoids do not generally show the erupting loops themselves but indicate 
the formation of the vertical current sheet below them that is the central element of the standard flare model. 
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1. Introduction 

Rising filaments, flares, and coronal mass ejections on the Sun 
often show the phenomenology of a loop-shaped magnetic flux 
system with fixed footpoints at the coronal base and signa- 
tures of magnetic twist. A single twisted magnetic flux rope 
appears to contain essential elements of the magnetic topology 
of the unstable, often erupting flux. Since the ratio of kinetic to 
magnetic pressure, the plasma beta, is very small in the inner 
corona,/? ~ lO"-' . . . 10"^, the flux system must be nearly force 
free in the quasi static pre-eruption state, i.e.. 




VxB = a(x)B . 



(1) 



Such major energy release events, which we will refer to 
collectively as flares in the following, often show a distinc- 
tive soft X-ray (SXR) brightening of sigmoidal (S or reverse-S) 



shape shortly prior or during their early, impulsive stages 


of de- 


velopment (e.g 


., Rust & Kumar. 1996, 


Marioharan et al.j 


199d 


Pevtsov et al. 




Il996l lAurass et all Il999l iMoore et all 


l200l[ 


Gibson et al.l 


2OO2I). The sismoidal shaoe is senerallv re 


jarded 



as a signature of coronal cuiTents in a twisted field and shows 
a hemispheric preference: forward (rev erse) S shapes dominat e 
in the southern (northern) hemisphere ('Rust&Kumail ll996l) . 
Typically the central part of such transient sigmoids is ap- 
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Fig. 1. Transient sigr noid observed iii an er uptive flare on 
1994 Oct 25 (from iManoharan etall 11994 in AR 7792 
which had primarilv lef t-handed twist jAurass et all 1 19991: 
Ivan Driel-Gesztelvi et aLi.i200Q) . 



proximately aligned with the neutral line of the normal pho- 
tospheric field component, often also with a filament. They 
are bright and more sharply defined in the middle, while their 
ends fan out more diff'usively at opposite sides of the neutral 
line. Sometimes their middle section is slightly split, giving the 
structure the appearance of two point-symmetric letters J. Each 
J may consist of several fibres (Fig.Q. 

The solar soft X-ray corona often also contains long-lived 
sigmoidal structures that appear to form a separate class. These 
typically exist for more than a day, do not show the concentra- 
tion in the middle, may be composed of many fibres that show 
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the S shape collectively rather than individually, and may form 
not only within but also between active regions. 

The sense of the sigmoidal shape and the handedness of the 
magnetic twist are correlated: IPevtsov et alJ ( fl997<) found that 
more than 90% of all sigmoids in those active regions that have 
a uniform sign of a in the photospheric magnetogram show 
forward (reverse) S shape for a > (<) 0, i.e., for right (left) 
handed twist. We will refer to this result as cr-a correlation, 
where cr > (<) stands for forward (reverse) S shape. The cor- 
relation is consistent with the hemispheric preference, since a 
is predominantly > (<) in the southe rn (northern) hemisphere 
(ISeehafer, 1990; Pevtsov et al., 1995). 

The transient sigmoids that brighten in flares (both confined 
and eruptive) have so far mostly been i nterpreted as an illumi- 
nation of the unstable flux system itself. iRust & Kur nar ( 1996) 
proposed that they show a sufficien tly twisted, ki nk-unstable 
flux rope. lPevtsov et alJ ( IT99 6). Moo re et ai] ( l200ll) . and others 
suggested further that the unstable flux rope acquired its super- 
critical twist through the reconnection of two nearly aligned, 
subcritical flux ropes near their end points, which conforms 
to the strong brightening in the middle of the sigmoids and 
may symbolically be writte n as J + J —> S . Transient sigmoids 
were thus regarded by Mo ore et alJ ( l200ll) as supporting their 
"tether cutting" model of flares (although they noted that an- 
other, much fainter structure was moving in sequences of SXR 
images in some of their events). 

The presence of a SXR sigmoid correlates with the 
likelihood of a flari ng active region to cause an eruption 
jCanfield et allfl999l) . A "standard" model for the main phas e 
of eruptive flares is now widely accepted (e.g., Shibata, 1999). 
It includes the formation of a large-scale vertical current sheet 
below the erupting flux, in which magnetic reconnection leads 
to the long-lasting energy release. This model is assumed to 
hold independently of the question whether the event is ini- 
tiated by an instability of a single flux rope, by reconnection 
between flux ropes (tether cutting), or by other mechanisms, 
e.g., "magnetic breakout" ( Antiochos et al. .199"^. 

The canonical instability of a twisted flux rope is the kink 
instability. In applications to the solar corona, the simplifying 
assumption of straight, cylindricallv symmetric flux ropes has 
always been used so far (e.g., lHoodlll992h . This is based on 
the typically large aspect ratio of solar magnetic loops. The 
instability occurs if the twist, O = lB^{r)/{rB,{r)), exceeds a 
critical value, Oq. Here I is the length of the flux rope, r is 
the (minor) radius, and B, and B^ are the axial and azimuthal 
field components, respectively. The critical value depends on 
the details of the considered equilibrium. 

In this paper we address the relationship between the kink 
instability of an arched twisted flux rope and the formation 
of current sheets and sigmoidal structures , starting from the 
anal ytical magnetic loop equilibrium by iTitov & DemouhJ 
(fl9'99). cited as T&D in the following (see their Fig. 2 for a 
schematic). This approximate, cylindrically symmetric, force- 
free equilibrium consists of a toroidal ring current of major 
radius R and minor radius a, whose outward-directed Lorentz 
self-force is balanced with the help of a field by two fictitious 
magnetic charges of opposite sign which are placed at the sym- 
metry axis of the torus at distances +L to the torus plane. That 



axis lies below the photospheric plane {z = 0} at a depth d. The 
resulting field outside the torus is current-free and contains a 
concentric magnetic X line between the torus and its centre. 
A toroidal field component created by a fictitious line current 
running along the axis of symmetry is added. This results in 
a force-free field inside the torus with the twist taking a finite 
value everywhere in the configuration, whose part in the vol- 
ume {z>0) is a model of a coronal loop (or flux rope). 

The presence of a toroidal field component turns the neigh- 
bourhood of the X line into a hyperbolic flux tube (HFT) which 
consists of two intersecting quasi-separatrix layers with ex- 
tremely diverging field lines (for a strict definition of HFTs, see 
iTitov e t al.U2002ll). The existence of the HFT is generic to such 
force- free loop configurations with a nonvanishing net current 
( Tito v et all Il999t) . which is important for understanding the 
origin of sigmoidal structures in twisted configurations. 

We present two characteristic cases of the kink instability, 
which show the different possibilities of current sheet and sig- 
moid formation in the T&D equilibrium, leaving a systematic 
study of the instability for a separate paper ( Torok et al., 20031) . 

T&D pointed out that the torus is also unstable with respect 
to global expansion (growing perturbations 6R > 0). This has 
recently been confirmed by Roussev et al. (2003). We expect 
that the global expansion instability leads to current sheets and 
sigmoidal structures in a manner largely analogous to the up- 
ward kink instability. 

We integrate the compressible ideal MHD equations using 
the simplifying assumption /3 = which has no influence on the 
qualitative evolution of the kink instability in its linear phase, 
during which the current sheets and sigmoids are formed: 



d,p = -V -ipu), 
p d,u - -p{u ■V)u + jxB , 
d,B = Vx(hxB), 
j = fi,'WxB. 



(2) 
(3) 
(4) 
(5) 



The equilibrium by T&D is used as the initial condition 
for the field, with the loop chosen to lie in the plane {x - 0}. 
The initial density distribution is specified such that the Alfven 
velocity is uniform, po = ^o/^o- Lengths, velocities, and times 
are normalized, respectively, by L, the initial Alfven velocity, 
v-dQ = Bo/(/U()Po)''^, and the Alfven time, = L/v-m- 

A modified Lax-Wendroff scheme is employed on a 
nonuniform Cartesian grid in a "half cube" [-L.v , L.v] x [0, Ly] x 
[0, Lj] with closed boundaries except the front boundary, {y - 
0), where line symmetry with respect to the z axis is imple- 
mented. The latter implies Ma,v(0, 0, z) = so that only vertical 
displacements of the loop apex are permitted. Viscosity is in- 
cluded in Eq. Q to impr ove numerica l stability. The numerical 
details are the same as in lTorok et alJ ( l2003i) . 



2. Kink instability 

The parameters of the basic loop equilibrium are chosen sim- 
ilarly to T&D: = L = 1 (50 Mm) and R = 2.2. The twist 
is specified via the number of turns that the field lines at the 
surface of the torus make about its axis (counted for the whole 
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Fig. 2. Top and side view of current density isosurfaces for <1) = 
4.9n. Left: \j\ - 0.15 y'max at f = 35; unperturbed case. Right: 
\j\ - 0.25 y'inax at f = 28; run with an upward initial velocity 
perturbation at the loop apex. |jc| < 1.5, |y| < 3, < z < 3. 



torus), taken here to be M = 15. This also fixes the minor ra- 
dius of the loop, a - 0.32. The resulting twist of the coronal 
part of the loop, averaged over it s cross se ction, 4> - A.9n, is 
about 1.40c of this configuration jTorok et al.. 2003 ) . Here the 
long- wavelength (m =1) kink mode develops spontaneously 
with downward displacement of the loop apex; it is initiated 
by the weak, downward-directed forces that result from the ini- 
tial discretization errors of the current density on the grid. The 
displacement and velocity of the apex rise exponentially un- 
til the apex height is reduced to k 2/3 of its initial value by 
t X 25, after which the instability starts to saturate. The con- 
figuration that results for left-handed twist {a < 0) is shown 
in Fig.|2(left) by current density isosurfaces. The kink pertur- 
bation forms a helical current sheet that is wrapped around the 
displaced parts of the loop at the interface with the surrounding 
plasma in a manner similar to the cylindrical case. The current 
density in this sheet rises exponentially as well, eventually ex- 
ceeding the current density in the loop. The projection of both 
the loop and the current sheet onto the bottom plane possess a 
reverse S shape, as is characteristic of sigmoids in left-handed 
twisted fields. If the handedness of the twist is reversed, then 
the m - 1 kink mode develops with identical time profile but 
reversed handedness of the loop perturbation, i.e., forward S 
shape of its photospheric projection. 

In order to study the kink mode with upward apex displace- 
ment, a small upward initial velocity perturbation is applied at 
the apex. It is ramped up over 5ra to - 0.01, then switched 
off. In this run the m - 1 kink mode develops with a similar 
growth rate as before, but again starts to saturate after t x 25, 
at an apex velocity of m, k 0.1. The resulting configuration in 
the case of left-handed twist is also shown in Fig. |2 (right). In 
addition to the helical current sheet that is again formed as the 
kink perturbation piles up the surrounding flux, a second cur- 
rent sheet, which has no analogue in the cylindrical case (but 
corresponds to the main element of the standard flare model), 
now occurs below the rising section of the loop. It is formed 
by a pinching of the HFT in the lateral stagnation flow set up 
below the rising loop apex. For a detailed analysis of the pinch- 
ing mechanism, which is a three-dimensional generalization of 
the well-known pinching of X points in the 2D case, we refer 



to lTitov etal] ( l2003l) and lGalsgaard et al.1 (l2003l) . The current 
density rises exponentially in both sheets, exceeding that in the 
kinked loop in the nonlinear stage of the instability. The rising 
unstable loop and the helical current sheet now both possess 
a forward S-shaped photospheric projection - opposite to the 
characteristic sigmoid shape in left-handed fields. Reversing 
the twist of the initial configuration leads to the upward m = 1 
kink mode with identical time profile but reversed handedness 
of the kink perturbation. It should be noted that the applied ve- 
locity perturbation is exactly vertical and so does not introduce 
a preference for the handedness of the kink perturbation. The 
two configurations of opposite sigmoidal shape shown in Fig.|2 
possess the same photospheric a. 

3. Formation of sigmoids 

In all four cases discussed in the previous section, the handed- 
ness of the kink perturbation of the loop axis (sign of "external 
writhe") equals the handedness of the field line twist (sign of 
"internal writhe"). This is a general property of the kink in- 
stability (kno wn as resona nce at the point of marginal stabil- 
ity; see, e.g. iLinton et all (l999j). Its immediate consequence 
is that transient sigmoids cannot be identified with rising kink- 
unstable loops. Since there are strong observational indications 
for the existence of substantial magnetic twist in eruptive flux 
systems, transient sigmoids must in general be associated with 
another structure in the unstable twisted flux. 

We follow the suggestion by T&D that bright SXR struc- 
tures outline magnetic structures with enhanced current den- 
sity, and hence enhanced dissipation, and calculate field lines 
that pass through the vertical current sheet below the upward- 
rising kinked loop (Figs.|3]and|4}. These field lines form a sig- 
moidal surface that not only follows the cr-a correlation but 
also corresponds nicely to further properties of transient sig- 
moids like the near-alignment with the photospheric neutral 
line, the concentration in the middle and the diffuse spread- 
ing at the ends. The occasionally observed splitting into two 
J-shaped parts, which may themselves consist of multiple J- 
shaped fibres, can readily be explained by assuming the exis- 
tence of multiple current density peaks within the current sheet. 
This can naturally be expected in a fibrous corona (we have in- 
dicated this effect in Figs. |3l and 0]by selecting field lines that 
pass closely to the z axis; in our simulation of a smooth ini- 
tial configuration, the current density peaks exactly at the z 
axis). Flux pileup in the stagnation flow at both sides of the 
current sheet may also lead to th e effect. Finally, cases of sig- 
moid lengthening, as described bv lManoharan et aD(ll996h . are 
consistent with the new interpretation only. 

Magnetic reconnection is expected to occur in the vertical 
current sheet, amplifying the brightening of the sigmoidal field 
lines at SXR and permitting transitions between continuous (S) 
and broken (double J) sigmoidal patterns (in both directions). 

With this interpretation, transient sigmoids cannot be 
considered as dire ct support of the tether cutting model 
dMoore et all 1200 ll) but, of course, they do not contradict the 
model in any manner. Our interpretation is not undermined by 
the fact that the kink-unstable loop does not succeed to evolve 
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Fig. 3. Field lines of the configuration shown in the right panels 
of Fig.|2l |x| < 0.8, lyl < 3, < z < 2.2. Loop field fines 
lie in tfie surface r - ajl (to separate tfie loop clearly from 
tfie sigmoids); one of tfiem is empfiasized to indicate tfie left 
fiandedness of tfie twist. Tfie remaining field lines pass tfirougfi 
two symmetrical vertical stripes close to tfie z axis tfiat bracket 
tfie pincfied HFT. Tfie normal component of tfie magnetic field 
in tfie bottom plane is sfiown color-coded (blue - positive, red 
- negative). Projections of tfie field lines are overlayed. 



a<0 




Fig. 4. Top view on selected field lines and on tfie surfaces tfiat 
tfiey form. |x| < 1, |y| < 3. Tfie field lines for a <Q are identical 
to tfiose plotted in Fig.|3l Tfie bottom panel sfiows similarly 
selected field lines of tfie system witfi reversed a at tfie same 
time. 



into a true eruption in our simulations dTorok et al.[l2003l) . be- 
cause transient sigmoids occur at tfie onset of eruptive events. 

Tfie downward-kinked loop may be associated witfi long- 
lived sigmoids. As tfie instability saturates, tfiis kink perturba- 
tion may lead to a stable structure wfiicfi follows tfie cr-a corre- 
lation and is nearly aligned in tfie middle witfi tfie pfiotospfieric 
neutral line. Tfie dip in tfie loop or in tfie flux below it supports 
tfie formation of a filament, wfiicfi is often observed in close 
spatial association witfi sigmoids ( Pevtsov. 2002). Furtfier stud- 
ies are required to cfieck tfiis conjecture. See iTorok & KliemI 
for an alternative interpretation of long-lived sigmoids 
as loops of enfianced current density but subcritical twist. 



4. Conclusions 

Tfie magnetic loop equilibrium by T&D is unstable against tfie 
long-wavelengtfi (m =1) ideal kink mode if tfie magnetic twist 
exceeds a critical value. Tfie fiandedness of tfie kink perturba- 
tion equals tfie fiandedness of tfie field line twist. Consequently, 
tfie pfiotospfieric projection of rising kink-unstable loops devel- 
ops a forward (reverse) S sfiape for left (rigfit) fianded twist. 
Tfiis is opposite to tfie strong observed {cr-a) correlation of 
sigmoid sfiape witfi magnetic fiandedness. Hence, transient 
SXR sigmoids, wfiicfi brigfiten in major energy release events, 
preferably in eruptive events, cannot in general sfiow tfie pre- 
sumed kink-unstable flux itself. We suggest tfiat transient sig- 
moids outline field lines wfiicfi pass tfirougfi tfie current density 
enfiancement in a vertical current sfieet formed below tfie rising 
loop. Tfiis interpretation leads to agreement witfi tfie cr-a corre- 
lation (and fiemispfieric preference) as well as witfi tfieir align- 
ment witfi tfie pfiotospfieric neutral line and witfi tfieir furtfier 
specific geometric properties. Transient sigmoids tfius indicate 
tfie formation of tfie vertical current sfieet tfiat is tfie central el- 
ement of tfie standard model for eruptive flares. Tfie question 
under wfiicfi conditions a rising kink-unstable loop develops 
into an eruption requires furtfier study. 

Downward-kinking loops take a sigmoidal sfiape tfiat fol- 
lows tfie observed correlation witfi fiandedness. Tfiey may be 
related to long-lived sigmoids and to filament formation. 
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